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ABSTRACT 

We present observations of the 22 GHz H 2 megamasers in the circumnuclear disk of NGC 3079 
obtained with the Green Bank Telescope. The data are analyzed for circular polarization due to the 
Zeeman-induced splitting of the H2O maser lines. No circular polarization is detected and we derive a 
la upper limit of 11 mG for the toroidal magnetic field at ~ 0.64 pc from the central black hole. This 
is the tightest upper limit for the magnetic field around a black hole to date. We use the magnetic 
field limit to derive an estimate of the mass accretion onto the central black hole. In addition to the 
polarimetric results, we present an observation of rapid variability in the maser lines, which we explain 
as weak interstellar scintillation. From the scintillation parameters, we estimate an intrinsic size of the 
mostly saturated maser features of ~ 12 microarcseconds. This is consistent with models assuming a 
thick, clumpy accretion disk. 

Subject headings: galaxies: individual (NGC 3079) — magnetic fields — masers-scattering 



1. INTRODUCTION 

Since AGN are thought to be powered by accretion 
onto a central super-massive black hole in which mag- 
netic fields likely play an important role (e.g. Blandford 
& Payne 1982), substantial magnetic fields are expected 
to be present throughout the entire accretion disk. Ob- 
■ servations of H2O maser circular polarization induced by 
Zeeman-splitting are an extremely useful method of prob- 
ing the magnetic field in regions where H 2 masers occur. 
Weak Zeeman splitting was initially observed in interstel- 
lar H2O masers using the Effelsberg 100 m telescope by 
Fiebig & Glisten (1989) and has since been used to de- 
termine the strength and structure of magnetic fields in 
, star-forming regions (e.g. Sarma et al. 2002; Vlemmings 
et al. 2006) and in the envelopes of late-type stars (e.g. 
Vlemmings et al. 2002). 

Recently, Herrnstein et al. (1998) and Modjaz et al. 
(2005) (hereafter M05) have used both the Very Large 
Array (VLA) and Green Bank Telescope (GBT) in an 
attempt to measure the circular polarization of the 
H2O megamasers in the circumnuclear disk of NGC 4258. 
The recent observations of M05 have placed stringent 
upper limits on the circumnuclear magnetic field of 
NGC4258. Here we present GBT observations of the 
H2O masers in the circumnuclear disk of NGC 3079 that 
were aimed at measuring the H2O maser Zeeman splitting 
due to the magnetic field in the disk. 

NGC 3079 is an almost edge-on spiral galaxy at approx- 
imately 16 Mpc (Sofue 1999). It shows either Seyfert-2 
or LINER activity and there is strong evidence from X- 
ray data that it contains an active galactic nucleus (AGN) 
(e.g. Irwin & Saikia 2003; Kondratko et al. 2005, and ref- 
erences therein). NGC 3079 also contains one of the most 
luminous 22 GHz H2O megamasers known to date (e.g. 
Henkel et al. 1984). Like that in NGC 4258, the H 2 maser 
emission of NGC 3079 has been interpreted as originating 
in a circumnuclear disk (e.g. Trotter et al. 1998; Sawada- 
Satoh et al. 2000; Yamauchi et al. 2004; Kondratko et 
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al. 2005) , although the exact disk models proposed dif- 
fer between authors. Kondratko et al. (2005) (hereafter 
K05) recently produced the first map of the full extent of 
the 22 GHz H2O maser emission of NGC 3079 (between 
Vlsr ~ 880 km s~ 1 and 1400 km s _1 ) and propose a model 
in which the accretion disk of NGC 3079 is thick, clumpy, 
flaring and undergoing star-formation. NGC 3079 is also 
the site of the first detected extragalactic submillimeter 
H2O masers, which also originate from the nuclear region 
(Humphreys et al. 2005). 

In addition to the measurement of Zeeman splitting, the 
GBT observations have allowed us to study the short-term 
variability of the masers. Previously, line flux variability 
of several tens of per cent on timescales of minutes has 
been observed for 22 GHz H2O megamasers of the Circi- 
nus galaxy (Greenhill et al. 1997; McCallum et al. 2005). 
As the timescale of intrinsic maser variations are at least 
an order of magnitude larger, the Circinus maser variabil- 
ity is likely due to a combination of weak and diffractive 
interstellar scintillation (McCallum et al. 2005). Here we 
present the first measurements of short timescale variabil- 
ity of the NGC 3079 H 2 megamasers and discuss the 
variability in the context of interstellar scintillation. 

2. OBSERVATIONS & DATA REDUCTION 

The observations of the H2O megamasers in NGC 3079 
were carried out at the NRAO 3 GBT between April 10 
and April 19 2006. At 22.2 GHz the full-width at half 
maximum (FWHM) beamwidth of the GBT is - 33", 
several orders of magnitude larger than the size of the 
H2O maser region, which is contained within approxi- 
mately 30 mas. The GBT spectrometer was used with 
a bandwidth of 200 MHz and 16,384 spectral channels. 
This resulted in a channel spacing of 0.164 km s _1 and 
a total velocity coverage of 2700 km s _1 , which was cen- 
tered on Vlsr = 1116 km s _1 . Furthermore, the data 
were taken with the dual-polarization receiver of the lower 
i\~-band (frequency range: 18.0-22.5 GHz) using the to- 
tal power nod observing mode. The two beams have a 
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fixed separation of 178.8" in the azimuth direction and a 
cycle time of 2 minutes was sufficient to correct for atmo- 
spheric variations. As a result, one beam of the telescope 
was always pointing at the source while the other beam 
was used for baseline correction. The total observing time 
was 30.2 hours of which 20.9 hours were spent on source. 
Once every <~ 1.2 hr point and focus observations were 
done on J0958+655. This source was also used as flux 
calibrator. The pointing corrections resulted in a pointing 
accuracy typically better than ~ 2". System temperatures 
ranged between 35 and 55 K, although an observing block 
of ~ 5 hr suffered from bad weather and as a result had a 
higher system temperature of ~ 65 K. 

For data reduction the GBT IDL 4 software package was 
used. The observations of J0958+655 yielded a gain-curve 
over the entire observation session which was linearly in- 
terpolated between calibrator observations and applied to 
the NGC 3079 data. Over the full observation run, the 
gain does not change by more than 10%, while the varia- 
tion over several hours in the flux density of J0958+655 is 
< 5%. Assuming that J0958+655 is not circularly polar- 
ized and that the circular polarization of the continuum 
emission of NGC 3079 is negligible, both the calibrator 
and the source data were used to determine a constant 
gain offset between the RCP and LCP data. It was found 
that the RCP gain needed to be multiplied by a factor 
0.945 to produce a Stokes V (= (RCP - LCP)/2) spec- 
trum without signal. To remove residual baseline varia- 
tions due to atmospheric and instrumental effects as well 
as removing the continuum emission from NGC 3079, we 
subtracted a baseline determined using the spectral chan- 
nels free of maser emission. We estimate the absolute flux 
calibration to be accurate to <~ 10%. As the magnetic field 
determination depends on the S/N and not the absolute 
flux levels, this uncertainty docs not influence the field 
strength determination. Finally, we produce source spec- 
tra in Stokes I and V, seen in Fig. 1, as well as RCP and 
LCP for the total combined data set. We find, from the 
spectral channels without signal, that we reach a la noise 
level of 1.4 mJy for I and V, while for RCP and LCP the 
noise level is 2.1 mJy. As the GBT observations presented 
here constitute hitherto the most sensitive observations of 
the H2O megamasers of NGC 3079, Fig. 2 shows a rescaled, 
hanning smoothed, total power spectrum highlighting the 
faint maser emission. The la noise level in this hanning 
smoothed spectrum is 0.9 mJy. 

3. RESULTS & DISCUSSION 

3.1. Maser Saturation 

Theoretical analysis in Nedoluha & Watson (e.g 1992) 
has shown that the degree of circular polarization of the 
22 GHz H2O masers depends on the level of maser sat- 
uration. The 22 GHz H 2 masers begin to saturate 
when R/T > 1 s" 1 and approach full saturation when 
R/T > 100 s _1 , where R is the rate of stimulated emis- 
sion and r is the maser decay rate. Following the analysis 
in M05, we find for the stimulated emission rate of the 
H 2 masers of NGC 3079 at 16 Mpc 

1 * - ' (1) 




Fig. 1. — Average total intensity (I, top) and circular polarization 
(V, bottom) spectrum of the NGC 3079 H2O megamasers. The 
solid vertical line at Vlsr = 1116 km s _1 indicates the systemic 
velocity. The dashed vertical lines indicate the maser features for 
which the variability parameters are given in Table 2. 
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Fig. 2. — Rescaled hanning smoothed total intensity spectrum of 
the NGC 3079 H2O megamasers highlighting the faintest detected 
maser flux. The solid vertical line indicates the systemic velocity. 
The dashed horizontal lines give the 3<r rms level of the hanning 
smoothed spectrum. 



Here we assume a cylindrical maser model with L the 
maser length and S v the maser flux density. The clumpy 
masing disk model from K05 predicts clump sizes be- 
tween 0.001 and 0.006 pc. Taking these sizes as lim- 
its on the maser length L, we find 70(5^/1 Jy) < R < 
2600(S'i,/l Jy) s _1 . The lifetime of the IR transitions in- 
volved in the masing process T = 1 s _1 , and the typical 
maser features are between 0.25 and 0.5 Jy; we thus con- 
clude that most of the H 2 masers of NGC 3079 are at 
least partially saturated and that several of the masers are 
even fully saturated. This is in contrast to the mostly un- 
saturated masers of NGC 4258 (M05) and to the prediction 
of R/T < 1 from numerical simulations by Watson & Wyld 
(2000). It is worth noting that the models from Watson & 
Wyld (2000) are specifically tailored for the thin accretion 
disk of NGC 4258 and are thus unlikely to be applicable 
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Fig. 3. — Average total intensity (I, top) and circular polarization 
(V, bottom) spectrum of the 925.5 km s" 1 blue-shifted H2O maser 
feature of NGC 3079. The dashed horizontal lines in the V-spectrum 
indicate the ±lcr errors. 



to the thick, clumpy accretion disk of NGC 3079. For high 
saturation, the spectral FWHM linewidth is expected to 
be Aw w 1.3 km s" 1 (e.g Nedoluha & Watson 1992; Vlcm- 
mings et al. 2002) for masing gas at intrinsic temperatures 
T ~ 300 K. In Table 2, we notice that although the ma- 
jority of the maser lines reach the expected width, several 
of the strongest lines have a Aw ~ 0.9 km s _1 . This either 
indicates that these lines are only partly saturated, that 
the intrinsic thermal line width of these lines is smaller, 
or a combination of these effects, thus highlighting various 
different conditions throughout the maser disk. 

3.2. The Magnetic Field 

As the H 2 maser lines in the spectrum of NGC 3079 
are very blended, we used the RCP-LCP cross-correlation 
method described in M05 to determine the Zecman split- 
ting and corresponding magnetic field strength for both 
the red- and blue-shifted masers. This cross-correlation 
method is specifically suited for blended maser lines of 
non-paramagnetic maser species, assuming the magnetic 
field is mostly constant across the spectrum, and it has 
been used to determine the upper limit of the magnetic 
field strength from the 22 GHz H2O maser polarization 
spectra of NGC 4258. As described in M05, the cross- 
correlation method yields the velocity shift q due to the 
Zecman effect, which is directly related to the magnetic 
field strength along the line of sight B^ by 

B \\ = ~i=^ • ( 2 ) 

Here Ap-F' is the coefficient that describes the relation 
between the circular polarization and the magnetic field 
for a transition between a high (F) and low (F 1 ) rota- 
tional energy level. The coefficient depends both on the 
intrinsic maser thermal width and the maser saturation 



level (Nedoluha & Watson 1992; Vlemmings et al. 2002). 
Here we have used Ap_pi = 0.020 km s _1 G _1 , as de- 
rived by Nedoluha & Watson (1992) for the merging of 
the three dominant 22 GHz H 2 maser hyperfine transi- 
tions. The results of this analysis are given in Table 1, 
where the rms error in velocity shift q, determined as de- 
scribed in the appendix of M05, yields the error in B\\. 
As the red-shifted masers are weak, the la error on the 
i?l|-field determination is large, and as no significant mag- 
netic field detection is made, we determine the 3a upper 
limit to be 807 mG. The error on the magnetic field in 
the blue-shifted masers is significantly smaller, but still 
no significant Zeeman splitting was detected. Here the So- 
upper limit is 33 mG, hitherto the tightest upper limit on 
the magnetic field in an accretion disk. 

In addition to the cross-correlation method for the heav- 
ily blended maser spectrum, an isolated H 2 maser line 
at Vlsr = 925.5 km s _1 also allows us to derive the mag- 
netic field using the regular Zeeman method (Vlemmings 
et al. 2002). Here the magnetic field is determined using 

P, Kiax ^min) „ . n> I \ ln\ 

v = ; = 2A F -f'B\\/Av 1 (3) 

J max 

where the fractional circular polarization Py is determined 
from the maximum and minimum circular polarization 
Knax and V min and the peak total intensity I max . The 
coefficient Ap-F 1 is described above and is again taken to 
be 0.020 km s" 1 G -1 , and Aw is the maser FWHM. As 
seen in Fig. 3, no significant detection is made in the cir- 
cular polarization spectrum. Using the calculation of the 
rms error limit as described in Vlemmings et al. (2006), we 
find a 3a upper limit of B\\ — 48 mG, similar to the up- 
per limit determined on the blue-shifted masers with the 
cross-correlation method. 

Assuming a standard accretion disk model, the magnetic 
field measurements at the blue and red-shifted maser fea- 
tures probe the toroidal magnetic field B$ at an average 
distance to the central black hole of ~ 0.64 pc. Since we 
found that the masers appear to be be either partially or 
fully saturated, the measured magnetic field strength could 
overestimate the true field strength by up to a factor of 4 
(Vlemmings et al. 2002, fig. 7). However, velocity gra- 
dients along the maser path can cause the true magnetic 
field strength to be underestimated by up to a factor of two 
(Vlemmings 2006). Additionally, the blending of maser 
lines in single-dish measurements can also lead to further 
underestimation of the actual field strength by a factor of 
two (Sarma et al. 2001). As we determined a similar upper 
limit on the isolated maser line at Vlsr = 925.5 km s _1 as 
on the blended blue-shifted masers, we adopt the conserva- 
tive 3a upper limit of B^ = 33 mG for the actual toroidal 
magnetic field strength. Using the model predictions for 
the magnetic field configuration in an accretion disk from 
Hawley et al. (1996) as described in M05, B^ dominates 
the total mag netic field and B tot w {W/IOY^B^. We thus 
take the total magnetic field strength to be Btot ^ 40 mG 
at a radius of 0.64 pc. Such a field strength value falls in 
the lower range of the magnetic field strengths observed 
for the H 2 masers in Galactic starforming regions (e.g. 
Fiebig & Giisten 1989; Sarma et al. 2002; Vlemmings et al. 
2006). Comparing to the limits found for NGC 4258 while 
taking into account the difference in the limit calculations, 
the magnetic field limit of NGC 3079 found in this paper 
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is ~ 2.5 times lower at 0.64 pc than that at 0.14 pc in 
NGC 4258. 

3.2.1. Mass Accretion Rate 

The mass accretion rate in the disk of NGC 3079 can 
be estimated from the magnetic field strength using the 
relation between the thermal pressure and the magnetic 
pressure within the accretion disk. Adapting formula 7 
from M05 to the case of NGC 3079 gives for the mass 
accretion rate M 

M< ID- 1 ' 5 f|f (fe) 2 (i^k) 1/2 

Here Tk is the temperature of the maser gas, r the dis- 
tance to the central black hole and Mbh the black hole 
mass. The ratio between gas and magnetic pressure is 
estimated from local three-dimensional MHD calculations 
to be (0) w 60 (Hawley et al. 1996). 

Thus for our total magnetic field upper limit of B tot = 
40 mG we find M < lO-^Mgyr" 1 at r = 0.64 pc, as- 
suming a smooth magnetic field configuration. As this 
estimate is based on the average magnetic field strength 
across the maser region, a much more complex magnetic 
field structure would lead to a higher upper limit by 
at most an order of magnitude. We can compare this 
upper limit to the estimate of the mass accretion rate 
found in K05 for their model of a thick, clumpy accretion 
disk. They find an accretion rate of 0.007Lboi,io ^yr 1 , 
where Z<boi,io is the AGN bolometric luminosity in units 
of 10 10 L©. With an estimate of the AGN bolometric lu- 
minosity from X-ray data of ~ 5 x 10 9 -5 x 10 10 L Q this 
indicates 10~ 25 < M < 10 _1 - 5 M Q yr _1 . Thus our re- 
sult for the mass accretion rate is fully consistent with the 
K05 model and implies an AGN bolometric luminosity of 
< 2.9 x 1O 1O L , which in turn implies an Eddington ratio 
of < 0.4. 

3.3. Variability 

Although the peak flux has varied by up to 50%, the 
overall shape of the total intensity spectrum in Fig. 1 is 
strikingly similar to the earlier spectrum from 2003 seen in 
figure 3 of K05. The only major difference is the appear- 
ance of a strong, somewhat isolated blue-shifted feature 
at Vjsr = 925.5 km s _1 . As shown in Fig. 4, this feature, 
in contrast with the other maser features, shows a steady 
decline in flux. Our time baseline is not sufficient to deter- 
mine a non-linear flux decrease. Using an error weighted, 
least squares fit to our individual observing scans of 4 min- 
utes yields a decrease in flux of 64 mJy/day. There are 
several possible mechanisms that can explain such gradual, 
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Fig. 4.— The light curve of the 925.503 km s" 1 blue-shifted 
H2O maser feature over the entire observation run. The thick solid 
line is a weighted least square fit to a simple linear decrease in flux 
density and corresponds to a decrease of 64 mjy per day. 

time- variable behavior (Maloney 2002). One possibility is 
a time- variable maser pump. The H 2 megamasers could 
be powered by X-ray emission from the AGN, as discussed 
by Neufeld et al. (1994). Variations in the AGN would 
then result in variations in the H2O maser flux. However, 
although the location of the Vlsr. = 925.5 km s _1 with 
respect to the other maser features is unknown, one would 
expect correlated flux variation in at least some of the 
other masers lines. Another more likely explanation is, 
that the variability is caused by the amplification of back- 
ground emission from another maser feature. As the disk 
of NGC 3079 is thought to be strongly inhomogeneous, 
the maser emission is narrowly beamed along the line-of- 
sight by the overlap of masing regions with similar Dopplcr 
velocities. As discussed in K05, chance alignments of mas- 
ing regions will then result in variability at timescales of 
severals days or longer. 

3.4. Scintillation 

In addition to the longterm variability seen in the fea- 
ture at Vlsr = 925.5 km s _1 , the maser features of 
NGC 3079 also vary on timescales of several tens of min- 
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a For the variability analysis the linear flux decrease has been removed. 
b Heavily blended feature. 





Fig. 5. — The variability of the blue-shifted H2O masers of NGC 3079. The spectra are averaged over 24 minutes, corresponding to 6 
observing scans. The average total intensity I spectrum from Fig. 1 has been subtracted for each of the spectra which are labeled with the 
MJD of the observation. Each spectrum is off-set by 0.5 Jy and the rms noise in the spectra is ~ 3.4 mjy. 



utes as seen in Fig. 5. Such rapid variability has been 
previously observed for the H 2 megamasers of Circinus 
(Greenhill et al. 1997; McCallum et al. 2005). As seen in 
Fig. 6, in which 7 hours of observation on the strongest 
red- and blue-shifted masers are shown, the observed vari- 
ations reach the level of 50% within the hour. This is 
similar to the flux variation seen in Circinus. Although 
small variations in time due to incorrect antenna gain cor- 
rections might still be present, the hour to hour gain vari- 
ations determined from the observations of J0958±55 are 
not expected to be > 5%. However, as the gain correc- 
tions using J0958+55 are only performed every 1.2 hours 
after correcting the telescope pointing, accumulated point- 
ing uncertainties and elevation effects can cause a global 
gain variation of up to « 10-25% between two consecutive 
observing blocks. While such large variations only occur 6 
times during our total observation run, the most extreme 
of these instances can be seen during the first hour of ob- 
servations presented in Fig. 6. These occasions are found 
to have no significant effect on our determination of the 



variability characteristics. An additional indication that 
the observed variability is not due to calibration errors is 
the fact that we find no significant correlation in the vari- 
ability of the different maser lines. 

To characterize the variability we have determined the 
discrete autocorrelation function (DACF) following Edel- 
son & Krolik (1988). The DACF of one of the maser 
lines is shown in Fig. 7. Following Rickett et al. (2002), 
we define the characteristic timescale for variability, T c h ar , 
as the time-lag where the DACF falls to 0.5. We fur- 
ther determine the depth of the first DACF minimum and 
the modulation index /i = <r/(S). Here (S) is the mean 
peak flux of the maser line and a the rms of the variabil- 
ity. The variability parameters are given in Table 2 along 
with the local standard of rest (LSR) velocity, Vlsr, and 
the full- width half- maximum (Aw). The velocity and Av 
are determined from fitting Gaussian profiles to the maser 
spectrum. However, many of the maser features are heav- 
ily blended, increasing the uncertainties in the velocity and 
Av determinations. We find that T c h ar ranges from ~ 900 



6 



Vlcmmings, Bignall & Diamond 



2.5 


: 




I 


925.503 km s 1 line - 


1.5 
1 


r 




• 

• 

• • 


• • - . _ 
• - • \ : 
. ••• "• * 






H 

.-■ 


I 1 h 


1 1 1 1 1 1 1 — - 

958.660 km s~ ! line - 


2.5 

g 2 


. *■**. 




••• 


'•. . • „/•-• = 

: 


1 1-5 












H 


I 1 h 


— i 1 1 1 1 1 1 — = 


* 14 








992.16 km s" 1 line ~ 


s 

E 1.2 






■ "•- 


.. . 

...»..•■ 

-•> 


0.25 




H 


I 1 h 


1 1 1 1 1 1 1 : 

1 190.65 km a" 1 line J 


0.2 
0.15 


[,8i 

- i 




*f 

I 


? i i i : 









2 


4 6 



t[hr] 



Fig. 6. — Light curves of the three strongest blue-shifted 
H2O masers lines of NGC 3079 as well as that of the strongest red- 
shifted line (bottom) observed on April 17th and 18th 2006. Data 
points arc plotted for every 4 minute observing scan. In the top two 
panels, the error bars are smaller than the symbol size. The typical 
la error bar is ~ 18 mjy. The light curve of the maser feature at 
925.503 km s _1 has been corrected for the linear flux decrease shown 
in Fig. 4. 
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Fig. 7. — Discrete, normalized autocorrelation function (DACF) 
of the 925.503 km s —1 maser line. The lags are binned in 4.5 min 
intervals. The DACF and the Icr error bars are calculated follow- 
ing Edclson & Krolik (1988). The solid vertical line indicates the 
characteristic variability timescale T^ar, where DACF= 0.5 and the 
dashed horizontal line is the first negative. 



to - 1600 S. 

As shown in McCallum et al. (2005), the variability of 



<~ 700 s of the Circinus H2O masers is likely due to a 
combination of weak and diffractive scintillation. How- 
ever, NGC 3079 is located at much higher galactic lat- 
itude (b = +48.36°) than Circinus, for which the line of 
sight passes through the galactic plane. From the Galactic 
electron density model NE2001 (Cordes & Lazio 2002), we 
find at the position of NGC 3079 that the predicted tran- 
sition frequency between weak and strong scintillation is 
vq = 8.8 GHz. The model uncertainties of NE2001 are sig- 
nificant at high latitude and for individual lines-of-sight, 
however NGC 3079 can reasonably be expected to be in 
the weak scintillation regime at 22.2 GHz. 

Using formula 6 from Walker (1998) for the frequency 
dependence of modulation index for a point source in the 
weak scintillation regime, we find from our observations, 
with a typical modulation index /1 = 0.1, the implied tran- 
sition frequency between weak and strong scintillation is 
vq = 4.4 GHz. The transition frequency would be higher 
if the maser size invalidates the point source approxima- 
tion. If the angular source size 9s is larger than the an- 
gular size of the first Fresnel zone 9 F at the distance of 
the scattering screen, then the modulation index is de- 
creased. At 22.2 GHz, 6 F = 59(0/lpc)" 1 / 2 ^tas, where z 
is the distance of the scattering screen. Assuming satu- 
ration, cylindrical masers are beamed by a factor of ~ 3 
(Elitzur 1994). As the clump sizes in the maser disk are 
estimated in K05 to be between 0.001 and 0.006 pc, this 
yields for the beamed masers, at 16 Mpc, an angular size 
between 4 and 25 /xas. As several maser features are only 
partly saturated, the beaming is more pronounced, and 
the angular sizes are even smaller. Thus, taking a typi- 
cal maser size of 63 = 10 /zas, the scattering screen would 
need to be at z ~ 150 pc to satisfy 6s/&f ~ 2.3 if the 
transition frequency is in fact close to the NE2001 model 
vo = 8.8 GHz. It is interesting to note that our calibration 
source J0958+655 (B0954+658), at an angular separation 
of 9.9° from NGC 3079, showed the prototypical extreme 
scattering event in 1980-81 (Fiedler et al. 1987), which was 
suggested by Fiedler et al. (1994) as being possibly asso- 
ciated with the edge of Galactic Loop III at a distance of 
145 pc. 

However, a distant screen implies a long scintillation 
timescale. At 22.2 GHz, the variability timescale for a 
point source is t F — 8100 V^z 1 / 2 , with Vism being the 
transverse speed of the scattering screen with respect to 
the source in km s^ 1 . For 6$ > Op, the characteristic 
variability timescale T c h ar = ^f^s/^f)- Thus for a screen 
at z = 150 pc and 6 S /6 F = 2.3, T char = 2.3 10 5 s" 1 . 
To reconcile this with the typical observed timescale of 
1200 s requires Vism ~ 190 km s _1 , which is much larger 
than the expected value of Vism ~ 40 km s _1 corresponding 
to the velocity of the Earth with respect to the LSR in the 
direction of NGC 3079 at the date of the observation. 

Alternatively, the observed rapid variations and lower 
transition frequency compared to the NE2001 model pre- 
diction could result from scattering in a nearby screen, as 
has been found for the handful of extremely rapid "intra- 
hour" scintillating quasars (Dennett-Thorpe & de Bruyn 
2000; Rickett et al. 2002; Bignall et al. 2003). Although 
there are some exceptions, the maser lines in Table 2 have 
a tendency to display the smallest T c h ar and largest \x for 
those lines with the lowest Av, which are thus thought to 
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be only partially saturated. This would be consistent with 
the corresponding maser features being the smallest due 
to the increased maser beaming. We thus expect 9s < Of 
for the partially saturated sources and expect 8s to be 
slightly larger than dp for the fully saturated sources. Tak- 
ing T c har ~ 1000 s from partially saturated masers to be 
ip, we find that the velocity and the screen distance have 
to satisfy z 1 / 2 = 0.12 Vism- For Vism « 40 km s^ 1 for 
a screen moving with the local standard of rest, this indi- 
cates z « 25 pc. This in turn gives 0s ~ &f = 12 /uas, thus 
the maser feature sizes expected from the scintillation are 
consistent with those estimated from the saturation level 
and the clump sizes from K05. 

Although there is significant estimation error in the 
DACF due to the limited sampling of the stochastic scin- 
tillation pattern, those masers with larger \i also tend to 
show deeper first minima in the DACF. This is also con- 
sistent with the sources with lowest fi, thought to be fully 
saturated, having 6s > Op- Rickett et al. (2002) showed 
that anisotropic scattering in a thin screen produces a deep 
first minimum or "negative overshoot" in the ACF, and as 
the source becomes more extended, this "negative over- 
shoot" is suppressed. A deep first minimum in the DACF 
is also seen for the intra-hour variable quasars and for the 
H2O masers in Circinus (McCallum et al. 2005), suggest- 
ing that anisotropic scattering in discrete "screens" is a 
widespread phenonemon. 

There is some indication of nearby scattering along other 
lines-of-sight close to NGC 3079 from the MASIV VLA 
Survey (Lovell et al. 2003). The MASIV Survey found 
that, while a large fraction of all compact, flat-spectrum 
extragalactic sources vary with modulation indices typ- 
ically in the range 0.01 to 0.1 at 5 GHz in a 72 hour 
period, only a tiny fraction (< 1%) have characteristic 
timescales of a few hours or less. J0949+5819, the clos- 
est MASIV source to NGC 3079 at an angular separation 
of 3.2°, showed variations with fi = 0.15 on a timescale 
of less than a few hours in 2002 January. This is likely 
to be due to scattering in a nearby screen within a few 
tens of pc. J0949+5819 was the most extreme variable 
observed in the first epoch of the MASIV Survey apart 
from the already known intra-hour variable J1819+3845 
(Dennett-Thorpe & de Bruyn 2000). Of the nine MA- 
SIV sources within 10° of NGC 3079, five sources showed 
significant variability in at least three out of four observed 
epochs (Lovell et al., in preparation), namely J0949+5819, 
J0946+5020, J0929+5013, J0958+4725, and our calibra- 
tion source J0958+655 which was also a previously known 
intra-day variable (IDV) (Wagner et al. 1993). While the 
fraction of nearby IDV sources is large but not exceptional, 
it is notable that the variations in J0929+5013 were also 
unusually rapid with a characteristic timescale less than 
a few hours. The presence of several rapid variables in 
the region of sky near NCG 3079 is suggestive of these 
sources being scattered by related structure in the very 
local Galactic interstellar medium. 

The proposed scintillation model could be tested by 
measuring the time delays and correlation between the 
maser scintillation patterns arriving at two widely sepa- 
rated telescopes. Such a two-station time delay experi- 
ment, for example between the GBT and Effelsberg tele- 
scopes, should also allow a determination of the anisotropy 
in the scintillation pattern and its position angle as well 



as the scintillation velocity (see e.g. Bignall et al. 2006). 

It is worth noting that short timescale variability of 
masers can also be caused by gain variation along the 
maser amplification path. As the characteristic timescale 
of such variability is ~ L/c, where L is the maser length, 
for the H2O maser emission of NGC 3079 this implies a 
variability timescale of this type of <~ 10 5 s. This is at 
least two orders of magnitude longer than the variability 
timescale detected here. 

4. SUMMARY 

The full velocity range of the H2O megamasers in the 
accretion disk of NGC 3079 has been observed using the 
GBT in both RCP and LCP. This has enabled us to de- 
termine a sensitive la upper limit of <~ 11 mG for the 
absolute poloidal magnetic field \B^\. This corresponds 
to a 3<7 upper limit for the total magnetic field in the ac- 
cretion disk of -Btot.iim = 40 mG at ~ 0.64 pc from the 
central black hole. This is the tightest magnetic field up- 
per limit in an AGN accretion disk to date. The magnetic 
field upper limit corresponds to a mass accretion rate of 
M < lO -1 - 7 M yr -1 , consistent with NGC 3079 accretion 
disk models and the X-ray luminosity. 

In addition to the magnetic field measurements, the 
H 2 maser lines of NGC 3079 were found to vary at char- 
acteristic timescales of the order of 1000 s. We conclude 
that this variability is due to weak scintillation caused 
by a scattering screen at ~ 25 pc. The scintillation 
timescales and modulation index imply maser feature sizes 
of ~ 12 /ias which corresponds to approximately 180 AU 
at the distance of NGC 3079. As the masers are found to 
be mostly saturated, maser beaming considerations imply 
that the masers occur in clumps with a size of <~ 0.002 pc, 
consistent with the model for a thick, clumpy accretion 
disk. 

This research was supported by a Marie Curie Intra- 
European fellowship within the 6th European Community 
Framework Program under contract number MEIF-CT- 
2005-010393. 
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